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Abstract 



< 

Two sharp existence and uniqueness theorems are presented for solutions of multiple vortices 
arising in a six-dimensional brane-world supersymmetric gauge field theory under the general 
gauge symmetry group G — £7(1) x SU(N) and with N Higgs scalar fields in the fundamental 
representation of G. Specifically, when the space of extra dimension is compact so that vortices 
are hosted in a 2-torus of volume |f2|, the existence of a unique multiple vortex solution repre- 
senting m , • • • , riN respectively prescribed vortices arising in the N species of the Higgs fields 
£ — . ' is established under the explicitly stated necessary and sufficient condition 
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£\i . where e and g are the J7(l) electromagnetic and SU(N) chromatic coupling constants, v mea- 

sures the energy scale of broken symmetry, and n — Y]j—y ni is the total vortex number; when 
the space of extra dimension is the full plane, the existence and uniqueness of an arbitrarily 
prescribed n-vortex solution of finite energy is always ensured. These vortices are governed 
by a system of nonlinear elliptic equations, which may be reformulated to allow a variational 
structure. Proofs of existence are then developed using the methods of calculus of variations. 

I Introduction 

The concept of solitons is important in quantum field theory. These static solutions, categorized 
into domain walls, vortices, monopoles, and instantons and often of topological origins [101 S3 HZ] , 
of gauge field equations, give rise to locally concentrated field configurations and are essential for 
the description of various fundamental interactions and phenomenologies. Vortices arise in two 
spatial dimensions and were first discovered by Abrikosov |T] in the form of a mixed state in a type- 

II superconductor in which the vortex-lines represent partial penetration of magnetic field into 
the superconductor as a consequence of partial destruction of superconductivity, in the context of 
the Ginzburg-Landau theory [21]. In quantum field theory, Nielsen and Olesen [43J showed that 



vortices arise in the Abelian Higgs model which may be used to model the so-called dual strings, 
better known as the Nambu-Goto |23] strings, which is a basic construct in string theory [63]. As 
a consequence of the Julia-Zee theorem [301 I54j . the vortex equations of the Abelian Higgs model 
in the static limit are exactly the Ginzburg-Landau equations. Although these equations are fun- 
damentally important, relatively thorough understanding of their solutions has only been achieved 
in a few extreme situations where one assumes either that magnetic field is absent [10^ [T4"l [371 E2] > 
that the solutions are radially symmetric [TB"! 144] , or that a critical coupling is maintained so 
that the interaction between vortices vanishes [29, 61J. In literature, this last situation is commonly 
referred to as the self-dual or BPS limit after the pioneering work of Bogomol'nyi [11] and Prasad- 
Sommerfield [45J. In fact, in the area of the study of non- Abelian gauge field equations, it is only in 
the BPS limit that tractable opportunities are available for gaining some fair understanding of the 
solutions of various equations, due to the difficulties associated with the presence of non-Abelian 
symmetry groups. The first array of existence results for non-Abelian vortices were obtained [551156] 
for the equations governing electroweak vortices formulated by Ambjorn and Olesen [21 S SI E], 
which were later sharpened [El [12]. Although the BPS vortices are present only when the coupling 
constants satisfy a certain critical condition (in the Ginzburg-Landau theory, this is the interface 
between type-I and type-II superconductivity; in the Abelian Higgs model, this is when the masses 
of the gauge and Higgs bosons coincide), such solutions exhibit a full range of elegant and unam- 
biguous features including exact topological characterization, quantization of flux and energy, and 
energy concentration, as anticipated from experimental facts. In recent years, the conceptual power 
of the BPS vortices in theoretical physics has been particularly witnessed in supersymmetric gauge 
field theory starting with the work of Seiberg and Witten [48] in an attempt to use non-Abelian 
color-charged monopoles and vortices to interpret quark confinement |38[ [39"1 |4"J~1 [571 [58l 160] . For 
surveys on this exciting topic, see [241 \3l\ |49"1 [52l 153] . Inspired by the importance of non-Abelian 
vortices in the linear confinement mechanism through a so-called dual Meissner effect, extending the 
classical Meissner effect in superconductors, some systematic research has recently been carried out 
aimed at understanding the various BPS vortex equations obtained in [151 \W\ [T7 1 125 1 1501 I5T1 152" ! [53] 
and a series of sharp existence theorems have been established [34 1 135 1 [36]. Unlike the BPS vortex 
equations in the electroweak theory [21 El [H [U El EH [56], the elegant structure of the BPS vortex 
equations in these supersymmetric models allow a complete understanding of their solutions, de- 
spite of the apparent sophistication associated with various underlying non-Abelian gauge groups. 
The main contribution of the present paper is to prove two sharp existence theorems for the BPS 
vortex equations arising in the supersymmetric U(l) x SU(N) gauge theory discovered in the work 
of Eto, Nitta, and Sakai [18] in the context of a six-dimensional brane-world scenario formalism 

[171 EH El] - 

An outline of the rest of the content of this article is as follows. In the next section, we 
review the non-Abelian multiple vortex equations of Eto, Nitta, and Sakai [18] and state our main 
existence theorems. In Section 3, we describe the nonlinear elliptic equation problem to be studied 
which is equivalent to the solution problem of the non-Abelian vortex equations. In Section 4, we 
consider multiple vortex solutions over a doubly periodic domain and identify a family of necessary 



conditions for the existence of such solutions. We then prove that these necessary conditions are 
also sufficient for existence. While doing this, we reveal various fine structures of the problem which 
will be seen to be useful for our later study of planar solutions. In Section 5, we prove the existence 
and uniqueness of a weak solution for the vortex equations over the full plane by the variational 
approach first developed for the scalar (Abelian) situation [29|. In Section 6, we first establish 
pointwise decay properties of solutions near infinity. We then strengthen these results by obtaining 
some exponential decay estimates. In Section 7, we use the obtained exponential decay properties 
of the solutions to establish various anticipated flux quantization formulas. 

2 Non-Abelian BPS vortex equations and main results 

In this section, we begin by a review of the non- Abelian vortex equations derived by Eto, Nitta, 
and Sakai in [18j . For convenience, we use x 1 and x 2 to denote the coordinates of extra dimensions 
in their six-dimensional gauge field theory. We will be sketchy and aimed at fixing notation since 
details can be found in [18] . 

We use {t a } to denote the generators of SU(N) (o = 1, • • • , N 2 — 1). Thus any element in the 
Lie algebra of SU(N), denoted by su(N), may be written as 

N 2 -l 



X= Y, xafa - (2- 1 ) 



o=l 



We use {qi} (i = 1, ••■ ,N) to denote N hypermultiplets in the fundamental representation of 
SU(N), i.e., each t& is a C -valued scalar Higgs field which is such that each of its value is taken 
to be a column vector. With (x £ ) = (x 1 ,^ 2 ), the gauge-covariant derivatives are given by 

V £ = d e -iA £ -iA e , £ = 1,2, (2.2) 

where Ag and Ag lie in the Lie algebras of U(l) and SU(N), respectively, i = y— 1, so that the 
induced gauge field strength tensors are defined by 

Fu., = d e A e , - dffAt, Fa, = d e A e , - d t .At - i[A e , A e ], £, £' = 1, 2, (2.3) 

where [•, •] denotes the matrix commutator. 

With the above notation, the non- Abelian vortex equations derived in [18] (see also [3 [26]) are 
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^ = yE^«' a = l,---,N 2 -l, (2.4) 

o / N 



V iqi = -V 2 q u i = l,---,N, (2.6) 



where e,g,v > are coupling parameters for which e represents Abelian (electromagnetic) gauge 
coupling, g non- Abelian (nuclear) gauge coupling, v the energy scale of the spontaneously broken 
ground state (vacuum), and t denotes the Hermitian conjugate. 



Now, following |18j . we collectively rewrite the Higgs fields and gauge fields in forms of N x N 
matrices, 

Q = (Qi,--- ,Qn), Ai = A e I N + A e , £ = 1,2, (2.7) 

where In denotes the N x N identity matrix. Then the ansatz 
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/N dmg{<p lr -- ,4> N }, (2.8) 

At = diag{-BJ,--- ,B?}, ^ = 1,2, (2.9) 

where fa are complex-valued scalar fields and B l = (Bj) are real-valued vector fields, i = 1, • • • , iV, 
further reduces the equations fj2.4|) — f |2T5T) into [18J: 

45L = 9 ^(v 2 - \fa\ 2 ) + (e 2 - £) ^ - ± £ |^| 2 ) , (2-10) 

(di+i^)^=i(Si+iS^, (2.11) 

for »= l,-- ,N. 

In view of [29J, we see from (|2.11|) that the zeros of each fa are isolated with integer multiplicities. 
We may use Z(fa) to denote the set of zeros of fa, 

Z(fa) = {Pi,i,--- ,Pi,nJ, i = l,---,N, (2.12) 

so that the repetitions among the points p^ s (s = 1, • • • , rii) take account of the multiplicities of 
these zeros. The first term on the right-hand side of (|2.10p indicates that these zeros enhance the 
'vorticity' field B\ 2 (i = 1, • • • ,N). However, the other terms on the right-hand side of (I2.10p 
complicate the situation so that we are not able to assert that the maxima of B\ 2 are achieved at 
Z(fa) (i = 1, • • • , N), which is what makes the problem interesting and challenging. 

There are two partial differential equation problems to be studied. 

The first one concerns the solutions of (|2.10p - (|2.1ip with prescribed zero sets given in (|2.12p 
over a bounded domain Q in M 2 so that the field configurations are periodic modulo the 't Hooft 
boundary condition |59t [6Tj at the boundary of fi. The physical relevance of this problem is that 
such solutions may arise when the extra dimensions are compactified as two-dimensional tori. For 
this problem, here is our result. 

Theorem 2.1. Consider the non-Abelian BPS vortex equations If2.10\) - i2.11\) for the field config- 
urations (fa,Bl) (i = 1, • • • , N ) over a doubly periodic domain £1 with the given sets of zeros stated 
in \2.12\) so that fa has exactly n« zeros in O (i = 1, • • • ,N). Then a solution exists if and only if 

the conditions 
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rij < \Q\ -\ 1 



n, i = !,-•• ,N, (2.13) 



are fulfilled simultaneously, where n = Xa=i n * * s ^ e total number of the N species of vortices. 
Furthermore, if a solution exists, it must be uniquely determined by the sets of zeros given in 
112.12]) , up to gauge transformations. 



Although the condition (|2. 13[) seems complicated, we may sum up the vortex numbers, rai, • • • , tin, 
to get a simple consequence of the inequalities stated in (|2.13[) , in the form 

2 2 

n<^-NM. (2.14) 

It is interesting to notice that, now, the non-Abelian coupling constant g does not enter the condition 
(I2,14p but the integer N which measures the 'size' of the non-Abelian symmetry. 

The second problem concerns the solutions of (|2,10p - (|2.1ip over the full plane M 2 . The form of 
the energy density described in [18] and requirement of finite energy naturally impose the boundary 
condition 

lim \<pi\(x) = v, i = l,-- ,N, (2.15) 

\x\— >oo 

for solutions. For this problem, here is our result. 

Theorem 2.2. The non-Abelian BPS vortex equations \2.10i) - l2.11\) over the full plane M 2 de- 
scribed by the field configurations ((pi, Bl) (i = 1, • ■ ■ , N) subject to the given sets of zeros stated in 
Ii2.12\) and the boundary condition \2.1Si) always have a unique solution. Moreover, such a solution 
realizes the boundary condition 112. 1 5\) exponentially fast. More precisely, for an arbitrarily small 
number e G (0, 1), there hold 

2 

\\^\-v\ + Y,\^- iB ^\ 2 + \ B ^\ =0 ^ {1 ' £) ^ lXl ^ i = h---,N, (2.16) 

£=1 
for \x\ sufficiently large. 

We note that, intuitively, Theorem 12.21 may be reinterpreted in view of Theorem 12.11 in the 
context of vortices in a domain of infinite volume. In such a situation, the condition (|2.13p of 
course becomes superfluous. 

Theorem 2.3. In both situations stated as in Theorems \2.1\ and \2.2\ the fluxes are quantized 
quantities given by the expressions 



I 



B[ 2 dx = 2irni, i = l,---,N, (2.17) 



which are seen to be determined by the algebraic numbers of zeros of the complex scalar fields (pi 
(i = 1, • • • , N ), with the integral evaluated over the doubly periodic domain Q or M. 2 , respectively. 

This result suggests that, although is not clear whether the vorticity field B\ 2 peaks at the zeros 
of the order parameter (pi, the total number of zeros of (pi determines the total vorticity (or flux) 
generated from the field B\ in the full domain. For this reason, we may still regard the locations 
of the zeros of (pi as the centers of vortices and refer to ni as the ith vortex number, which extends 
the concept of vortices in the classical Abelian Higgs model |29j . As a consequence, the integer n 
defined in (I2.13P is well justified to be called the total vortex number of the solution. 



With the quantized fluxes given in (|2.17p . the total non-Abelian vortex (string) tension, Tna> 
may be computed [18] to assume the elegant exact value 



N 



T NA = ^E/^ d - = ^E^ ( 2 - 18 ) 

The above theorems will be established in the subsequent sections. 

3 System of nonlinear elliptic equations 

To proceed, we now adapt the complexified variables and derivatives defined by 

z = x 1 + ix 2 , B i = B i 1 + iB i 2 , d = ^(d l -id 2 ), d = ^(d 1 + id 2 ), (3.1) 

and convert the equations (|2.1U|) - (|2.1ip into the following system 

4i(<W - m) = ^(|^| 2 - , 2 ) + (e 2 - £) f 1 £ |^| 2 - vA , (3.2) 

BhKffi = l -B\ (3.3) 

away from the possible zeros of <f>i, for i = 1, • • • ,N. 

Inserting (|3.3p or B l = — 2i<91n<^j into (|3.2p (i = 1, • • • , iV) and using the relation A = <9 2 + <9| = 
439 = 4:dd, we arrive at the equations 



Alnl&l 2 



^(ltf-^)^(^4)(l||tf-^), (3.4) 



away from the zeros of <j>i, i = 1, • • • ,N. Thus, with the notation in (|2.12p and the new variables 

Ui = bx\<f>i\ 2 , i = l,---,N, (3.5) 

we obtain the following system of nonlinear elliptic equations 



1 / 2\ / I N \ n i 

^-v 2 ) + ^[e 2 - 9 -ji-^ U '-v 2 J+^Y, 6 P^ x ^ * = 1 r--,N, (3.6) 



governing the interaction of N species of vortices located at the prescribed set of points 

Z = uf =1 Z(4> i ), (3.7) 

which are the set of zeros of the complex scalar fields </>i, • • • , (J>n- 

We shall consider the solutions of the system (|3.6p over a doubly periodic domain (a 2-torus) 
Q, realized by the 't Hooft periodic boundary condition |59|. IBT] and over the full plane R 2 . In the 
latter situation, we need to observe the boundary condition (J2.15P at infinity. That is, 

limui = 21nt>, i = l,---,N. (3.8) 

| re I — >oo 

In the next section, we shall first concentrate on the doubly periodic situation. 



4 Necessary and sufficient condition for doubly periodic solutions 

In this section, we study the equations (|3.6|) defined over a doubly periodic domain Q. We conve- 
niently rewrite these equations as 



N 



Aui = Y^ Oij (e nj - v 2 ) + 4vr ^ S Pie (x), i = 1, ■ ■ ■ , N, 



8=1 



where 





a ij ~ 


1 

JV 


(f 


9 2 

2N 


Let 


u® be a solution to 










An? 


= - 


4irrii 

\n\ 


+ 4t 


• [6j 


) The substitutions 









+ 5,, 



13 2AT' 



i,i = !,••• ,7V. 



n i 



s=l 



i = l,... ,iV. 



Ui = Ui +Ui, i = !,■•• ,JV, 



recast the equation (|4.ip into 



JV 






-''-« a ) + ^ i = l,-,*. 



We use boldfaced letters to denote column vectors in W 1 . Thus, we set 
U=(U 1 ,--- ,U N y, G= (e M ? +c/l ,--- ,e< +u »y, 



Ifil 






l^l 



u 2 ^ ojvj 



and let A = (a 



ijJNxN 



be the N x N matrix defined by (14. 2ft . Then 



,4 



/ a + b a a 
a a + b a 



\ 



a a 



a \ 
a 

a + bj 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



where a = j? ( \ — j^j and b = ijjj which should not be confused with the group indices used in 
Section 2. 

Now the equations (14. 4p can be written in the vector form 



AU = ,4G + F. (4.5) 

This system looks difficult to approach. In order to tackle it, we shall seek for a variational principle. 



To find a variational principle, we need to use the property of the matrix A. It is easy to check 
that the matrix A is positive definite and its eigenvalues are 



\ 1 = Na + b= —, X 2 



Xm = b = — . 

2N 



(4.6) 



Then, by the Cholesky decomposition theorem |22j . we know that there is a unique upper triangular 
N x N matrix T = (tij) for which all the diagonal entries are positive, i.e., tu > 0, i = 1, • • • , N, 
such that 



A = T T T. 



(4.7) 



In fact, by direct computation, we have 
in = Va + b, tu = ti 3 = ■ ■ ■ -- 



hN = — = "1 > 0, 
til 



*22 = v/(a + 6) -Of, t 23 =t 2A 



t 2 N 



a — on 



t-22 



a 2 > 0, 



t 



(N-1)(N-1) 



\ 



N-2 



(a + b)-J2 a l * 



N-2 o 



(N-l)N 



i=l 



EiV— z 
i=l a 

£(AT-l)(iV-l) 



= ajv-l > 0, 



tviv 



N 



7V-1 



(a + 6)-^af. 



i=l 

-TT"\ — 1 



Set v = (t>i, • • • ,vn) t , L = (T r ) = (lij)isfxN- We introduce the new variable vector 

v = LU or U = L _1 v = T T v. 



Then (14.511 takes the form 



Av = TG + LF. 



(4i 



(4.9) 



With the convention oln = 0, we may write (|4.9p in the component form 

JV i 

Avi = t u e u ° +tiiVi+ ^tJi «*«* + a . J^ e u ° +t ^+^i=i a * Vk + ^ Ujfj, i = l,--- ,N. (4.10) 

j=i+l jr' = l 



It is easy to check that the above system of equations (|4.1Up are the Euler-Lagrange equations 
of the functional 



Setting 



J ( v ) = / U E i v ^i 2 + E ^ +i ^ +E - a ^ + E E w,- 

^ [ i=l i=l i=l W=l 

% = f e u ° +tiiVi+ ^i=^ akVk dx, i = 1, • ■ • , N, 
Jo, 



>dx. (4.11) 



(4.12) 



(note that these g^'s should not be confused with the N scalar "quark" fields in the fundamental 
representation of the gauge group denoted by the same notation in Section 2) we may integrate 
flUDD to obtain 

N i 

tuqi + <Xi } ] qj = —\£l\ ^ kjfj = Pi, * = 1, ■■■ ,N. (4.13) 

j=i+l j=l 

Since ta > (i = 1, • • • , N) and en > (i = 1, • • • , N — 1), the definition of qi (i = 1, ■ • ■ , N) given 
in (|4.12p and the relation (J4.13P lead to the necessary condition 



Pi >0, i = !,-•• ,iV, 



(4.14) 



which appears to be complicated. 

In order to arrive at an explicit form of the necessary condition, recall the structure of the 
matrix T given in (|4.7p . Thus, with p = (pi, • • • ,Pn) t and q = (<Zi, • • ■ , </tv) t , we can rewrite (J4.13P 
in the matrix form 

Tq=-\n\LF = p. (4.15) 

With L = (T T )"\ we can solve (I4TT5D to get 



\n\T- 1 LF = — |n|A _1 F. 



(4.16) 



On the other hand, since for any invertible N x N matrix D and the column vectors X and Y 
in M, N satisfying Y r D -1 X ^ 1, the matrix M = D — XY T is invertible and 



M- 1 = (D- XY T )" 1 = (I + [1 - Y T D- 1 X\- 1 D- 1 XY T )D- 1 . 
Applying the formula (J4.17P to the matrix 

A = diag{6, • • • , 6} - (-a, • • • , -a) r (l, • • • ,1), 
we have 



A' 1 



1 



/(7V-l)a + 6 



b(Na + b) 



-a (N-l)a + b 



\ 



Inserting (I4.19P into (14.16J) . we obtain the solution 



—a 



N 



-a \ 



(N-l)a + bJ 



(4.17) 



(4.18) 



(4.19) 



where 



|fi| I b fi b(Na + b)^ fj 

2|„, 47TO ATT . 

AT 

n = y^n^ 



(4.20) 



(4.21) 
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is the total vortex number. Substituting the values of the constants a and b in terms of the coupling 
constants e and g in (|4.20p . we have 

«i = v 2 \n\+8v^- 1 ^n-^n i >0, i = l,--,N, (4.22) 

which lead us to the necessity of the condition (12. 13ft . 

Below, we shall show that, under the condition (|2.13|) . the system (|4.10p has a solution. We 
will use a direct minimization method as in [34J. Furthermore, in order to gain more insight to 
the technical structure of the problem, we shall also sketch a constrained minimization method to 
approach the problem. 

We use 14 7l,2 (Q) to denote the usual Sobolev space of scalar- valued or vector- valued O-periodic 
L 2 -functions whose derivatives are also in L 2 (Q). In the scalar case, we may decompose W 1,2 (Q) 
into W l > 2 (Vl) = R © W l ' 2 {Q) so that any / G W l ' 2 (Vt) can be expressed as 

/ = / + /, /et, f€W^(n), [ fdx = 0. (4.23) 

It is useful to recall the Moser-Trudinger inequality [BJ [19] 

J e u dx< Cexpf-^- f |Vu| 2 dx) , «6^' 2 (fl). (4.24) 

With (|4.24|) . it is clear that the functional I defined by (14. lip is a C 1 -functional with respect to its 
argument (vi, • • • ,Vn) £ W ' 2 (Q), which is strictly convex and lower semi-continuous in terms of 
the weak topology of W 1 ' 2 ^). 

We can suppress the functional / given in (14.1ip into the form 

„ / N N \ N 

J(v) = / l~J2 |V^i| 2 +Je^ + S=i a "* 1 dx - ^PiVi. (4.25) 

•^ n V i=l i=\ / i=l 

Applying the Jensen inequality and using the fact that £j, > (i = 1, • • • , iV) and a, > 
(i = 1, • • • , N — 1), we have 

J e "°+*»fe+**)+U=Vfe & +^)dx > (njexp f_ _L J u o dx \ eX p ( ^ + ^ « fe ^ j 

= ffi e*«2* + nii Q ^, » = 1, . . . , TV. (4.26) 

Substituting fCMD into (ET251) . we have 

7 ( v ) " o / Y2 \ Vi, i\ 2 dx > V aie^^+S^i Qfe ^ - V M r (4.27) 

2 /o — ^ — J — ^ 

Ji ' i=l i=l i=l 

In order to control the lower bound of ()4.27p . we recall the relation Tq = p with the quantities 
qi (i = 1, • • • , N) as defined in (J4T22D . That is, 

AT 

Pi = tiiQi + Oi ^2 Qj, i = 1, • ■ ■ , N. (4.28) 

j=i+i 
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Therefore we have 



Now set 



Then we arrive at 



N N / i-\ \ 

53jhv» = E *( ^ + X] ° fc ^ • ( 4 - 29 ) 

i=l 8=1 V fe=l / 



8-1 
fc=l 



tei + ^ajfcUfc, * = !,••• ,iV. (4.30) 



N N 



J ( v ) - \ I E i v ^i 2 dx ^ E (^ - M£) • ( 4 - 31 ) 

"^ 8=1 8=1 



Thus, using the elementary inequality 

< ce 6x - ax, a, 6, c> 0, i£K, (4.32) 

in (|4.3ip . we obtain the coercive lower bound 



£ f i- ta 



JV JV 



I(v) - £ /" X>*| 2 ds > f> (l + ln N) 



(4.33) 



It follows from (|4.33p that if (v) is bounded from below and we may consider the following direct 
minimization problem 

77 = inf{J(v)| ve W 1,2 (fi)}. (4.34) 

Let {(vi , • • • , v\r )} be a minimizing sequence of H4.34n . Since the function 

F(u) = ae u - qu, (4.35) 

where a, q > are constants, enjoys the property that F(u) — > oo as u — > ±oo, we see from (|4.3ip 
that the sequences {wf 1 } (i = 1, • • • , JV) are all bounded where w\ n is defined by (|4.30p by setting 
W-i = Uu an d 2Z« = V4 {i = I; ' ' ' > JV, n = 1, 2, • • • ). Inverting the transformation (I4.30p . we see 
that the sequences {v^ } (i = 1, • • • , JV) are also bounded. Without loss of generality, we may 
assume 

v^ 1 — > some point v 4 £l as n — > 00, z = 1, • • • , JV. (4.36) 

In addition, using (|4.33[> . we conclude that {Vv\ } (i = 1, • • • , JV) are all bounded in L 2 (f2). 
Therefore, it follows from the Poincare inequality that the sequences {■?;■ } (i = l,--- , JV) are 
bounded in VF 1,2 (f2). Without loss of generality, we may assume 

v\ n —7- some element v^ G VF ' (fi) weakly as n — )• 00, i = 1, • • • , JV. (4.37) 

Obviously, ^ oo) G W 1>2 (n) (t = 1, • • • , JV). Set t;J oo) = t), W +ij°° ) (i = 1, • • • , JV). Then (jQ5j| 
and (|4.37p lead us to see that f 2 — >■ t^ weakly in PF ' 2 (Q) as n — > 00 (i = 1, • • • , JV). The 
weakly lower semi-continuity of I enables to conclude that (v^ , • • ■ ,v^° ) solves (|4.34p . which is 



12 



a critical point of I and a classical solution to the system (|4.10p in view of the standard elliptic 
theory. 

Since the matrix A is positive definite, it is easy to check that the functional I is strictly convex 
in W 1,2 (Q). So it has at most one critical point in W 1,2 (ft), which establishes the uniqueness of 
the solution to the equations (14.101) . 

Below, we will develop our methods further by presenting a constrained minimization approach 
to the problem. For convenience, we rewrite the constraints (14.121) collectively as 



J 4 (v) = I e u ° + '* m+ £fc=i ° kVk dx = qi , i = 1, • • • , N. 
Jn 



(4.38) 



Recall that the values of q\ , • • • , q^ are given by (14.22P which are obtained by solving the system 
of equations f)4. 13|) , or 

N i 

tm + on Y, Qj + M Yl Wi = °> * = X > " " " ' N - ( 4 - 39 ) 

j=i+l j = l 

We consider the multi-constrained minimization problem 

n = inf{/(v)| v G W 1 ' 2 ^), Jx(v) = <?i, ■ ■ ■ , Jjv(v) = q N }. (4.40) 

Suppose that ()4.40p allows a solution, say v = («!,••• , fjv). Then there are numbers (the 
Lagrange multipliers) in R, say £i, • • • , £jv, such that, for i = 1, • • • , N, 



f Ivvi- Vvm + tiie u ° +tiiVi+ ^i^ a ^ +a t Y e u °i +t ^ + ^ k => kVk + J^hjfj 

r - N r 

= &u / e^+^+^i^u-idx + a, V / e^ +t ^ + ^ k ~^ akVk Widx, (4.41) 

where w±,--- ,wjy are test functions. Letting w\ = ■ ■ ■ = wn = 1 in the above equations and 
applying f|4.39j) . we arrive at 

N 

j=i+l 

Consequently, £i = • • • = £at = so that (|4.4ip is exactly the weak form of the system (|4.10p . 
In other words, the Lagrange multipliers disappear automatically and a solution of (|4.40p solves 
(|4.10p . Hence, it suffices to find a solution to (I4.40|) . 

In order to approach (|4.40p . we use the notation (J4.23P to rewrite the constraints (J4.38P as 

e^fii+Efci «kv k f e «?+t«i>i+Etei a kv k dx = q h i = 1, • • • , N, (4.43) 

Jn 

which may be resolved to yield 

i 

v, = ^ Zij (In ^- -lnJ,(v)), i = l,---,N, (4.44) 

3=1 
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where v = (vi,- ■ ■ ,vn) t and lij are the entries of the lower triangular matrix L = (T T ) with 
T = (Uj). 

To proceed, we use the constraints (I4.38D to rewrite the action functional (14. lip or (I4.25P as 



1 N . N N 

i=l Jn i=l i=l 

N i N / i \ 

= ^ J^ p^ - In J,- (v) + J^ I % - k Y^ h ln 9j > ( 4 - 45 ) 

i=i i=i i=i \ j=i / 

where we have inserted (|4.44[) . However, using the relation q = L T p, we can rewrite (J4.45P as 

N N 

J ( v )-2^/ |V^| 2 dz = J>lnJi(v)-0), ( 446 ) 

t=l "^ i=l 

where Co is a constant depending only on L, p, and q. By virtue of the Jensen inequality, we have 

Ji(v) > |fi|exp( / u°dx) =o-i, i = l,---,iV. (4.47) 

Using the condition q± , ■ ■ ■ , qx > and the lower bound (|4.47[) in (|4.46p , we obtain the following 
coercive inequality 

1 N r N 

/(v)--^ / |Vw 4 | 2 da;> J^Wi-Co- (4.48) 

1=1 ^ 8=1 

Now the proof of solvability of (|4.40p follows from a standard argument. 

In fact, let {(v\ , ••• ,v N )} be a minimizing sequence of (|4.40|) . In view of (J4.48J) and the 
Poincare inequality, we see that {{vi , • • • , v^ )} is bounded in W 1,2 (fi). Without loss of generality, 
we may assume that {(v\ , • • • , ujy )} converges weakly in W 1,2 (fJ) to an element («i, • • • , Ujv). The 
compact embedding 

^ ll2 (fi) -»• L p {n), p > 1, (4.49) 

implies that (?)},••• , vjy ) — >• (i?i,--- , wjv) in L p (il) (p > 1) as n — > oo. In particular, v_ t = 
(i = 1, • • • , 2V). In view of (|4.24|) and (|4.49|) . we see that the functionals defined by the right-hand 
sides of (|4.44p are continuous in fy (i = 1, • • • ,N) with respect to the weak topology of W ' 2 (Q). 
Therefore, vf 1 — > some v^ £ R (i = 1, • • • , iV) as n — )■ oo, where u, is given in (|4.44j) . In other 
words, v = («!,••• , Vjv) = fei + vi,--- ,Utv + vn) satisfies the constraints (J4.38P and solves the 
constrained minimization problem (|4.4U|) . 

5 Solution on full plane 

In this section, we prove the existence and uniqueness of the solution to equations (|4.ip - (|4.2p over 
M? satisfying the natural boundary condition m = 21nu (i = 1, • • • ,N) at infinity as given in (|3.8[) . 
Under the translation m >— > Ui + 2 ln v (i = 1, • • • ,N) and the rescaling of the coefficient matrix, 



14 



v 2 ciij i-> dij (i,j = !,-•• ,N), the equations (j4.ip become 



N 



A Ul = y, M e% - l ) + 4vr Yl <w (»). » = i, • • • > w, 



where 



8=1 



v ( i 9 \ P 5 w ■ • 
o« = — — e — — -f Oii — rr; »i 7 = 1» " • , iv, 

J 2JV V N J 2N 



subject to the boundary condition 

Ui — > 0, i = 1, • • • , iV, as |x| — >• oo. 
As in |29|, l62j 163] , we introduce the background function 



Ui(x) = -^ln(l + /u|x-p iiS | 2 ) , /i>0, i = 1, 



,iV. 



s=l 



Then we have 



/ij(x) = 4 E s =i (At+ | a; _p iiS p)^ 
i = i,--,jv. 

Let Ui = u® + Ui, i = 1, ■ ■ ■ ,N. Then the equations (I5.ip become 



A? 



Set 



AC/, = ^ a 4J (e^+^ " 1) + hi(x), i = 1, • • ■ , N, 

3=1 



U = ([/!,-•• ,U N ) T , G = (e< +Ul -1,--- , e u °N +UN -1) T , 
H = (fci(x),..- ,Mx)) T . 



Thus, the equations (|5.6[) can be written in the vector form 

AU = AG + H, 

where A = (a,ij) NxN . 

As before we use the transformation (j4.8|) to change (|5.7j) into 



i«i+*ii«<+SJ, = l "fc^fc 



JV 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



Avi = tu \e u l +tiiVi+ ^=i akVk - 1J + on Y, f e tt § + *" 1, ' +s *=i a ^ fc - l) + J^ hjhj, (5.8) 

j=«+i j=i 

for i = 1, • • • , JV. It is direct to check that (|5.8|) are the variational equations of the energy functional 

N 



J(v) = fit) ^ Vv ^ + (e"° +4 "^+^=i^^ - e u ? 



t-l 
fc=i 



+ ftfj > dx, (5.9) 
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where 

i 

gi = J2kjhj, i = l,---,N, (5.10) 

i=i 
which may also be rewritten as 



i-1 

tuVi + ^2 a k v k 
fe=l 



J(v) = Y^ \ -llVwilH + e u se f ""' + S=i at " fc - 1 

i=l IV L /c=l J / 2 

+ |e u ?-l,^i + ^a fc x; fe J + (gi,Vi)A, v 1: - ■ ■ ,v N e W 1 ' 2 ^ 2 ), (5.11) 

where (•, -)2 and || • ||2 denote the inner product and norm of L 2 (M. 2 ), respectively. 

It is clear that the functional I is a C 1 -functional with respect to v and its Frechet derivative 
satisfies 

DI(y)(v) = f; I IIV^IH + (t llVl + ^ t!)fcl e"? +1 "»' + Si« - 1 J + ( 9i , Vi ) 2 1 . (5.12) 

i=l I V fc=l / 2 > 

Set 

i-1 

Wi = tnVi + y^akVk, i = l,--,N, (5.13) 

fc=i 

w = (wi, • • • , wn) t , and g = (gi, • • • ,qn) t - Then w = T T v or v = Lw. Thus 

N „ „ N 

y\{9i,Vi)2= g r vdx= / h T L T Lwdx = S2(H i ,w l ) 2 , (5.14) 

^ J ®> 7k2 ^ 

where H = (Hi, ■ ■ ■ ,H N ) T = L T Lh = A^h. Inserting (J57T3J) and ([5TT4"]) into (pU2j) . we obtain 

AT 

I>/(v)(v) = J^ {||V^||i + (^,e<^ - 1 + #,) } . (5.15) 

«=i 

To estimate the right-hand side of (|5.15p . we consider the quantity 

M(w) = (w,e U0+w -1 + H) 2 , (5.16) 

where w, uo, H stand for one of the functions Wi, u®, Hi, for i = 1, ■ ■ ■ ,N, respectively. 

As in [29], we decompose w into its positive and negative parts, w = w + — W- with w + = 
max{w,0} andw_ = —mm{w,0}. Then M(w) = M(w+) + M(w-). Using the inequality e 4 — 1 > t 
(t 6 R), we have e uo+w — 1 > w + uq, which leads to 

M(w+)> \\w+\\ 2 2 + (w+, u + H) 2 > ^\\w + \\l--\\u + H\\l, (5.17) 



where we have used the fact that H,uq £ L 



2(m2\ 
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On the other hand, using the inequality 1 — e * > t/(l + 1) (t > 0), we can estimate M{—W- 
from below as follows: 

M(-ii/_) = (w-,l - H - e uo + e uo [l - e- w -}) 2 



> [w^,l-H-e u ° + 

„2 



W- 



wz 



1 + W- 



-,1-H + 



W- 



-,«<) 



l + w_ 



-,1-H-e 



«o 



(5.18) 



From the definition of the function i?, we may choose /i > large enough so that iJ < 1/2 
(say). Note also that H, 1 - e u ° E L 2 (M 2 ). Thus, we have 



iu_ 



l + w_ 



,1-H-e' 



»o 



< 



»r 



1 + W- 



■dx+ \\1-H -e 



«0||2 



(5.19) 



Summarizing these facts, we see that (|5.18[) enjoys the lower bound 



M(-w-) > 



1 



wz 



l + W- 



dx-C, 



(5.20) 



where and in the sequel, C denotes a generic but irrelevant positive constant. 
From (J5TT1) and (fOUj) . we have 



M(w) > 



Using (|5.2ip in (|5.15p . we arrive at 



1 



w 



1 + ku 



dx - C. 



(5.21) 



^ i N (■ 2 

™(v)(v,-£llv»,ll!>i£/ M ,T^-c 



(5.22) 



Moreover, since the matrix T is invertible and v and w are related through w = T r v, we can 
find a positive constant Cq such that 



A' 



JV 



^llV^|||>CoX)ll V ^lli- 



(5.23) 



i=l 



i=l 



Inserting f)5.23|) into (|5.22j) . we get 



£/(v)(v)>CoEll V ^ll2 + 7E/ TrK dx - C - 
We now recall the standard Gagliardo-Nirenberg-Sobolev inequality [32j [33] 
/ f 4 dx<2J fdx [ \Vf\ 2 dx, feW 1 ' 2 ^ 2 ). 



(5.24) 



(5.25) 
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Consequently, we have 
fdx^j = 

< 2 

< 4 



l/l 



l + l/l 

7 2 

(i + l/l) 2 

7 2 
(i + l/l) 2 



(1 + 1/1)1/1 dx 

dx/ (/ 2 + / 4 )dx 

JR 2 



da; 



/ 2 dx 1+ / |V/| 2 dx 



1 
< - 

~ 2 



J 2 dx + C 1 + 



f- 



i'i 



(l + l/l) 1 



■ dx 



+ 



n4 N 



|V/| 2 dx 



As a result of ()5.26j) . we have 



/ 2 dx 



<C 1 + 



|V/| 2 dx + 



/ 2 



dx 



(5.26) 



(5.27) 



(l + l/l) 2 

From (|5.24p . (|5.27p . and the relation between v and w, we may conclude with the coercive lower 
bound 



/ N N \ 

DI(y)(y) > d I £ Nh + Yl II V ^H 2 



C 2 , vi, 



,v N 6W rl ' 2 (R 2 ) 



(5.28) 



where C±, C<i are some constants. In view of the estimate (|5.28|) . the existence of a critical point of 
the functional / in the space M /1,2 (IR 2 ) follows in a standard way. 
In fact, from (|5.28j) . we may choose R > large enough such that 

inf{l>I(v)(v)|v=(vi,... ,^)e^(i 2 ), \\v\\ w i,2 {m = R} > 1. (5.29) 

Consider the optimization problem 

r] = inf {l(v)| ||v|| w i,2 (M 2) < R} . (5.30) 

Let {v' 71 -*} be a minimizing sequence of (|5.30|) . Without loss of generality, we may assume that 
this sequence is also weakly convergent. Let v be its weak limit. Thus, using the fact that the 
functional / is weakly lower semi-continuous, we have I(v) < ij. Of course ||v||^-i,2m2) < R because 
norm is also weakly lower semi-continuous. Hence I(y) = r\ and v solves (|5.30p . We show next that 
v is a critical point of the functional /. In fact, we only need to show that v is an interior point, 
or ||v||jyi,2(]R2) < R. For suppose otherwise that ||v||ypi,2(R2) = R. Then, in view of (I5.29h . we have 



lim 
t-+0 



/(v - tv) - J(v) 



d 



/(v - tv) 



t=0 



(DI(y)) (v) < -1. 



(5.31) 



t dt~ 

Therefore, when t > is sufficiently small, we see by virtue of (15.311) that I(v — tv) < /(v) = r\. 
However, since ||v — tv||^i,2m2) = (1 — i)R < R, we arrive at a contradiction to the definition of v 



or (15.30p . Thus v is a critical point of i". 

Finally, the strict convexity of I says that I can only have at most one critical point, so we have 
the conclusion that / has exactly one critical point in W 1,2 (M. 2 ). Of course, this critical point is a 
solution of (I5.ip . which must be smooth by virtue of the standard elliptic regularity theory. 

The asymptotic behavior of the solution will be studied in the next section. 
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6 Asymptotic behavior of planar solution 

Let v = («i, • • • , vn) denote the solution of (|5.ip obtained in the previous section. Here we aim to 
establish the pointwise decay properties for v. Our tools are based on elliptic L p -estimates and the 
maximum principle. 

To proceed, we first recall the following embedding inequality [20l [29j [33] 



LP(R2) < I IT 



, s P-2 

L2 



" f V 1 > 2 (K 2 )> P > 2 - (6- 1 ) 



From (|6.ip . we may infer e^-l£ L 2 (IR 2 ) for / G W 1,2 (M. 2 ). To see this fact, we use a MacLaurin 
series expansion to get 



N 2 ^ 2 s - 2 
e J 

' , a' 

s=3 



s-2 



l ) =f 2 + E-^f s - ( 6 - 2 ) 

By virtue of (|6.ip and (|6.2p . we obtain 

i^ - i|ii 2(K2) < \\f\\h m + E ^ (t 1 ^) 2 il/il^w < 6 - 3 ) 

which confirms our claim that e^-l£ L 2 (M 2 ) because it is easily shown that the series on the 
right-hand side of (16. 3p is convergent. 

Now consider (|5.8j) . Since Uj G W 1,2 (M 2 ) (z = 1, ••• , AQ and right-hand side of (|5.8p may be 
rewritten as 

iv 
t .. e ^ Ltim+Z&ccM _ i) + £.. ( e «? - l) + ^ E { e "° (e^+Eti *"* - l) + (e u ° - l) } , 

i = l,---,iV, (6.4) 



we see that the right-hand side of each of the equations in (|5.8p belongs to L 2 (R 2 ). Thus we 
may resort to the standard elliptic L 2 -estimates to deduce that V{ G W 2,2 (E?) (i = 1, • • • , AT"). In 
particular, v(x) — )• as |x| — > oo because we are in two dimensions. 

We can establish similar decay properties for \Vvi\ (i = 1, • • ■ , N). To this end, we rewrite the 
right-hand sides of (|5.8p as 



All these belong to the space L P (R ) for any p > 2 due to the embedding W 

and the definition of uf. Therefore elliptic L p -estimates enable us to conclude that V{ G W 2,P (M?) 

(i = 1, • • • ,N;p>2). In particular, we have \Vi>i\ — > as \x\ — > oo, i = 1, • • • , AT. 

To obtain suitable exponential decay estimates for the solution, it suffices to consider (I5.ip 
outside the disk Dr = {x G M 2 | \x\ < R} where 

R > max{|pi s | \i = 1, 2, • • • , N, s = 1, 2, • • • , m}. (6.6) 



tu L^ _ A e Um+j:i = \ a k v k + t .. ^t.-^+EUl ^ k v k _ {• 

+ai Y, { ( e " J ° " l) e'»^+^*=i Qfc,;fc + te"' + £*=i Qfc ^ - l) } , » = 1, • • • , AT. (6.5) 
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For convenience, we write (|5.ip in Mr \ Dr in the form 

N N 

Au i = ^2a ij u j + ^2a ij (e u J -Uj -1) , i = 1,2, • • • ,N. (6.7) 

3=1 3=1 

With the vectors u = (u\, • • • , un) t and V = (e Ul — u% — 1, • • • , e UN — u^ — 1) T , we may rewrite 
(I677D as 

Au = ,4u + AV. (6.8) 

Let O be an N x N orthogonal matrix so that 

T AO = diag{Ai,A 2 ,--- ,X N }= A, (6.9) 

where Ai, A2, • • • , A^ are as given in (|4.6|) . Thus, in terms of the new variable vector 

U = (C/i,--- ,U N ) T = O t u, (6.10) 

Substituting (|6.10p into (|6.7p and using (|6.9p and the behavior of U — > as \x\ — > 00, we arrive at 

AU 2 > 2AjvU 2 - a(z)U 2 , i£l 2 \ D R , (6.11) 

where a(x) — > as |x| — > 00. Consequently, for any e 6 (0, 1), we can find a suitably large R £ > R 
such that 

r2 n o\„ 11 itt2 ™^-Tn>2 



AU 2 > 2A^v (1 - -J U 2 , x g M 2 \ D^. (6.12) 

Thus, using a suitable comparison function, the property U = at infinity, and the maximum 
principle, we can obtain a constant C(e) > to achieve 

u 2 = U 2 (x) < C(e)e- (1 - e)v ^l x l, |x| > i? £ . (6.13) 

We next derive some exponential decay estimates for |Vttj| (i = 1, • • • , N). For given £ = 1,2, 
we differentiate (|6.7p to obtain 

N 

A(d eUi ) =Y / a ij {x)e u '{d t u j ), i = l,---,N. (6.14) 

3=1 

Set v = (dgui, ■ ■ ■ ,d(UN) T and E(x) = diag{e U1 ^', • • • ,e UN ^ x '}. Then the system (|6.14|) becomes 

Av = Av + A(E(x) - I N )v, (6.15) 

where In is the N x N unit matrix. Consequently, we have 

Av 2 > 2v r Av 

= 2v r ^v + 2v T A(E(x) - I N )\ 

= 2 (O r v) T A (O r v) + 2v T A{E(x) - I N )v 

> 2X N v 2 - a(x)v 2 , x£R 2 \D R , (6.16) 

where a(x) — > as \x\ — > 00. Hence, as before, we obtain the estimate 

N 

^2\d e u i \ 2 = v 2 <C(e)e-^- £ '> V2X ^^, \x\>R £ , 1 = 1,2. (6.17) 
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7 Consequences of asymptotic estimates 

For the function hi given in (|5.5p . we can directly compute to get 

hi dx = A-nrii , i = l,---,N. (7.1) 



Furthermore, for the solution (U%, • • • , Un) of the system (|5.6p . we have U{ = Ui — U® where u® 
is defined by (J5.4J1 (i = 1, • • • , JV) and (u±, • • • , un) is the unique solution of (|5.ip subject to the 
boundary condition (J5.3J) whose derivatives are seen to vanish at infinity exponentially fast. Thus 
we infer that |V£/i|,--- , |VC/jv| all vanish at infinity at least as fast as |x| -3 . Consequently, we 
have 

/ AUidx = 0, i = I,--- ,N. (7.2) 

Jr 2 
Thus, integrating (|5.6p and applying (|7.ip and (|7.2p , we obtain the quantized integrals 

N 

y^an (l-e^)dx = 4vrn,, i = l,---,N. (7.3) 

Pi ^ 

Now recall the relation between the functions u\, • • ■ , un and the Higgs scalar fields <j>\, ■ ■ ■ , 4>n 
described in Section 2. Since 

|&| 2 = v 2 e*S i = l,---,N, (7.4) 

we see that 

|^| 2 - v 2 = v 2 (e Ui - 1) = o(e-^-^^^^), (7.5) 

if \x\ is large. Using (|7.5p in (|2.10p . we conclude that the curvatures B\ 2 (* = lj""" >^0 vanish at 
infinity exponentially fast at the same rate. Besides, since the Higgs fields <f>i,--- ,4>n and gauge 
fields Bj, ■ ■ ■ , B¥ may be constructed from u\, • ■ ■ , u^ following the expressions 



4>i{z) = vex.pl -Uj(z) + i^arg(x 

V 3=1 



(z-Pi, s )\, (7.6) 

B\{z) =r -Re{21d In (j)i(z)}, B\(z) = -Im{2i51n<^(z)}, (7.7) 

z = 1, • • • , N, we see that the covariant derivatives obey 

2 2 

J2 K d t - Mt)<t>i\ 2 = ^\Vu t \ 2 = o(e-( 1 - £ )^ 2 ^l a; l), i = l,---,N, (7.8) 



when \x\ is sufficiently large. Hence, inserting the value Ajv = 5 2 /2-/V, we obtain all the decay 
estimates stated in Theorem 12.21 

Finally, applying (|7.3p and (|7.4p in (|2.10p and noting (|5.2p . we obtain the quantized flux formulas 
stated in Theorem 12.31 in the full plane case. 

In the situation of a doubly periodic domain, the same flux quantization conclusion follows 
simply from integrating (|4,ip and no further consideration is necessary. 
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